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The charmonium yields are expected to be considerably suppressed if a deconfined medium is 
formed in high-energy heavy-ion collisions. In addition, the bottomonium states, with the possible 
exception of the T(1S) state, are also expected to be suppressed in heavy-ion collisions. However, in 
proton-nucleus collisions the quarkonium production cross sections, even those of the T(15'), are also 
suppressed. These "cold nuclear matter" effects need to be accounted for before signals of the high 
density QCD medium can be identified in the measurements made in nucleus-nucleus collisions. We 
identify two cold nuclear matter effects important for midrapidity quarkonium production: "nuclear 
absorption", typically characterized as a final-state effect on the produced quarkonium state and 
shadowing, the modification of the parton densities in nuclei relative to the nucleon, an initial- 
state effect. We characterize these effects and study the energy, rapidity, and impact-parameter 
dependence of initial-state shadowing in this paper. 



I. BASELINE TOTAL CROSS SECTIONS 



To better understand quarkonium suppression, it is necessary to have a good estimate of the expected yields. How- 
ever, there are still a number of unknovifns about quarkonium production in the primary nucleon-nucleon interactions. 
In this section, we discuss models of quarkonium production and give predictions for the yields in a number of collision 
systems. 
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TABLE I: For each ion species at the LHC, we give the maximum beam energy per nucleon and the corresponding beam 
rapidity. Using the maximum proton or deuteron beam energy: Ep — 7 TeV and j/p = 9.61; -Ed = 3.5 TeV and yd ~ 8.92 
respectively, we present the maximum center-of-mass energy per nucleon; rapidity difference, y^^ = yt — Va (i = p, d); and 
center-of-mass rapidity shift, AyJ^ ~ vIub/'^j for P + ^) d+A and A + A collisions. Note that there is no rapidity shift in the 
symmetric A + A case. 



Since the LHC can collide either symmetric {A + A) or asymmetric {A + B) systems, we present results for p + p, 
p + A, d+A and A + A collisions. We consider d+A collisions since the d-|-^ center-of-mass energy is closer to the 
v4 -|- v4 collision energy than top energy p + A collisions. The maximum ion beam energy per nucleon is the proton 
beam energy, Ep = 7 TeV, times the charge-to-mass ratio, Z/A, of the ion beam. Thus the maximum deuteron beam 
energy is half that of the proton beam, = 3.5 TeV. The ion beam energies are given on the left-hand side of TableU 
for five reference nuclei: oxygen, g^O; argon, fgAr; krypton, IgKr; tin, ^J^Sn; and lead, ia^Pb. Note that we use the 
average elemental A since a sample may contain an admixture of several isotopes of different A. 

In addition to the A -f A center-of-mass energy, we also show the maximum p + A and d+A per nucleon center-of- 
mass energies, ^/s'^ = yj/AEp^ (iEa- Because -Ep,d is typically greater than Ea, the center-of-mass rapidity can shift 
away from y — 0. The total shift is y^^g = Ui — Ua (* = P: d) while the center of mass shifts by half this amount, 
Ayl^ = TableUshows the maximum nucleon-nucleon center-of-mass energy per nucleon, the rapidity difference 

between the two beams, y^^, and the center-of-mass shifts ior p + A and d+A collisions. (The Z/A ratio is the same 
for d and O thus Ay'^'-' = 0.) Only ^s„„ is given for symmetric A + A collisions since there is no rapidity shift. 

If there were no cold nuclear matter effects on the production cross sections at a given energy, the per nucleon 
cross sections would all be equal. However, the nuclear parton distributions (nPDFs) are known to be modified with 
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respect to the free proton PDFs as a function of parton momentum fraction x. At low x, x < 0.05 (shadowing region), 
and high x, x > 0.2 (EMC region), the nuclear structure function, F^{x), the weighted sum of the charged parton 
distributions, is suppressed relative to that of the deuteron, i^^(x), while, in the intermediate x region, the ratio 
2F^/AF2 is enhanced (antishadowing) in nuclear deep-inelastic scattering (nDIS). We refer to the modification of 
the parton densities in the nucleus as a function of A, x and /x^ in general as shadowing. While a combination of nDIS 
and Drell-Yan data can separate the nuclear valence and sea quark densities, there is no direct probe of the nuclear 
gluon density, rather it is inferred from the scaling violation. 

Gluon fusion dominates quarkonium production up to xf ^ 0.7 already at fixed-target energies [H, including over 
the entire accessible rapidity range at the LHC, see Fig.[T] Thus while the modification of the gluon distributions in 
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FIG. 1: (Color online) For M = 4 and yTsj;;; = 20 (solid), 40 (dashed), 200 (dot-dashed), 5500 (dotted) and 14000 (dot-dot- 
dot-dashed) GeV, we give the average value of the center-of-mass rapidity, y, in p + p collisions as a function of xp- Note the 
absolute values on y and xf: in the center-of-mass frame, \xf\ < 1 and \y\ < j/max so that the curves shown here are diagonally 
reflected around xp ~ y ~ 0. 



nuclei is the most important for quarkonium studies, it is unfortunately the most poorly known. There are, however, 
a number of indirect constraints on the gluon density. The scale evolution of and momentum conservation provide 
two important constraints. Most of the low- a: nDIS data are at relatively low scales, below the minimum scale of a 
number of PDF sets and therefore less useful for studies of perturbative evolution. RHIC data on hadron production 
are an exception since intermediate pt hadron production occurs at relatively low x and at perturbative scales. At 
relatively high x, the shape of the PHENIX midrapidity n° data 0] helps pin down the nuclear gluon density in the 
EMC region. 

Quarkonium production occurs at sufficiently large scales to provide further constraints on the nuclear gluon PDFs. 
There are some drawbacks however: the quarkonium production mechanism is not fully understood, even in p + p 
collisions, and the energy dependence of nuclear absorption is not well known. In the remainder of this section, 
we discuss the quarkonium yields in various collision systems; the implementation of modified PDFs for the nuclear 
parton densities; and quarkonium absorption by nucleons. 

Early studies of high energy quarkonium production, particularly at high pT, were performed in the context of the 
color singlet model (CSM) which calculates direct production of a quarkonium state with definite total spin, parity and 
charge conjugation. The CSM predicted that the Xci state, produced directly from gg fusion, would have a much larger 
cross section than direct color-singlet J/tp production which requires a 3-gluon vertex Instead, measurements of 
direct J/^p and Xc production showed that the J/-0 cross section was, in fact, larger than the Xc cross section [3]. 
However, the CSM can describe charmonium production in cleaner environments such as photoproduction Q and 
more recent formulations of the CSM Q can reproduce the magnitude of the J/-0 data at RHIC Q. The different 
kinematics of the modified CSM probes larger values of x and thus reduces the shadowing effect 

Nonrelativistic QCD (NRQCD) is an effective field theory in which short-distance partonic interactions produce QQ 
pairs in color singlet or color octet states which then evolve into a quarkonium state, as characterized by nonperturba- 
tive matrix elements @. The first term in the NRQCD expansion is equivalent to the CSM. The octet contributions 
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are sufficient to explain the J/tj} yield at the Tevatron. However, the NRQCD approach has so far failed to describe 
quarkonium polarization . 

Perhaps the simplest approach to quarkonium production is the color evaporation model (GEM) which treats heavy 
flavor and quarkonium production on an equal footing. The quarkonium production cross section is some fraction, 
Fc, of all QQ pairs below the HH threshold where H is the lowest mass heavy-flavor hadron. Thus the CEM cross 
section is simply the QQ production cross section with a cut on the pair mass but without any contraints on the 
color or spin of the flnal state. The color of the octet QQ state is 'evaporated' through an unspecified process which 
does not change the momentum. The additional energy needed to produce heavy-flavored hadrons when the partonic 
center of mass energy, -s/j, is less than 2mH, the HH threshold energy, is nonperturbatively obtained from the color 
field in the interaction region. Thus the quarkonium yield may be only a small fraction of the total QQ cross section 
below 2mH- At leading order, the production cross section of quarkonium state C in an A + B collision is 

"^"^^^2^^ ^ FcY, C^" I dx^dx2 [ dz'dz 
d\d^b T^J^^^l J J 

X ff{xi,p^,r,z) ff{x2,fJ.'^,b-f,z') aij{s) S{s - xiX2S^„) , (1) 

where A and B can be any hadron or nucleus, ij = qq or gg and aij{s) is the ij — >■ QQ subprocess cross section. If one 
or both of the collision partners, A and B, is a proton, then the transverse, r, and longitudinal, z, spatial parameters 
may be replaced by delta functions, J d'^rdz6{r)6{z), and the parton densities are simply //^(xi, , r, z) = ff-{xi, p?). 
Our calculations use the NLO QQ code of Mangano et al. jll| with the 2raH mass cut in Eq. ([!]) and use the same 
parameters as in Refs. [H, with the MRST parton densities optimized to obtain agreement with the QQ 
cross section, as described in Ref. [l^. The factor Fc can also be fit with other parton densities such as CTEQ6M 
[l6j . When the same mass and scale parameters are used, the energy dependence of the cross section is very similar, 
see e.g. Ref. [isj . 

All these formulations: CSM; NRQCD; and CEM assume the validity of collinear factorization which relies on the 
separation of the initial and final states. Collinear factorization was proven to be effective at all orders for the Drell-Yan 



process some time ago [17] . A subsequent paper by Collins, Soper and Sterman showed that the factorization process 
was correct for heavy fiavor production up to corrections of order (1/M) [l^. Thus while higher-order corrections to 
the charm cross section are large, collinear factorization is generally assumed to hold and, indeed, the scale dependence 
of the approximate NNLO-NNLL charm cross section is seen to stabilize and the next-order corrections are not as 
large [19]. Higher- twist effects that might signify factorization breaking, such as intrinsic charm, are generally most 
important at forward rapidities in the light-cone formulation [20l . [2l| . Factorization has not been strictly proven for 
quarkonium where the final quarkonium state may be connected to the initial state by soft gluons. In the CSM, the 
color singlet matrix element is derived from quarkonium decays where the initial state plays no role. Factorization is 
most difficult to prove for NRQCD. It depends on the universality of the nonperturbative matrix elements. However, 
recent works have shown that a redefinition of these matrix elements allows factorization to be restored [1^ [2^ . The 
CEM is closest in spirit to the calculation of the open heavy flavor cross section so that collinear factorization should 
work equally well in the two approaches. Since our calculations are in the CEM, we use collinear factorization to 
calculate quarkonium production at the LHC. 

To go beyond p + p collisions, the proton parton densities must be replaced by those of the nucleus. Then the 
collision geometry and the spatial dependence of the shadowing parameterization also need to be considered. We 
assume that if A is a nucleus, the nuclear parton densities, ff^{xi,is^,f,z), factorize into the nucleon density in 
the nucleus, pA{r,z), independent of the kinematics; the nucleon parton density, /f(a;i,/i^), independent of A; and 
a shadowing ratio, Sp c^{A,xi, p.'^,r, z) that parameterizes the modifications of the nucleon parton densities in the 
nucleus. The first subscript, P, refers to the choice of shadowing parameterization, while the second, S, refers to the 
spatial dependence. Thus, 

f.f{xi,fi'^,r,z) = pA{s)Sps{A,xi,p,'^,f,z)ff{xi,fi'^) , (2) 
ff{x2,p^,b-f,z') = pB{s')S^ps{B,X2,fi'^,b-r,z')fJ{x2,fi'^) , (3) 



where s = Vf^ + and s' = y \b — rp + z'^. 

The nucleon densities of the heavy nucleus are assumed to be Woods-Saxon distributions ^23] and are normalized so 
that / d^rdzpA{s) = A. With no nuclear modifications, g(A, x, Q'^,r,z) = 1 and integration of the nuclear parton 
densities over the spatial variables gives 

d^dhdzdz' /^(xi,/x2, r, z)f^{x2,pM- r, z') = ABff{xi,p^)fP{x2,fi^) . (4) 
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The impact-parameter averaged shadowing parameterization measured in nDIS is recovered by integrating Sp^s over 
the volume, weighted by the nuclear density, 



d'^rdzpA{s)S'p^s{A,x,^l^,r,z) = Si,iA,x,^l^) . (5) 



We discuss more details of the spatial dependence of S'p,s in Section 2.3. Most available shadowing parameterizations, 
including the ones used here, ignore the small effects in deuterium. However, we take the proton and neutron numbers 
of both nuclei into account. The impact-parameter integrated up and down quark distributions, needed for the qq 
contribution to quarkonium production, are calculated as 

f^{x,^,') = {ZAS|,^{A,x,^l^)fp{x,^,^) + NAS'p^iA,x,^,^)f:;{x,^l^)) (6) 

for q = u and d, assuming that, as for the proton and neutron parton densities, S'p„ — Sp^ and 5'p„ = Sp^ and 
similarly for the antiquarks. 

To obtain the rapidity distribution from the total cross section, an additional delta function, d{y — 0.5 ln(a;i/a;2)), is 
included in Eq. ([T]). At leading order, the parton momentum fractions xi and X2 are simply xi.2 = (y^s/s^^ ) exp(±j/). 
In this notation then, in the forward rapidity region oi a p + A collision, a;i, the proton momentum fraction, is larger 
and X2, the parton momentum fraction in the nucleus, is smaller than the midrapidity value, x — ^s/s^^ . 

Some of the uncertainties in the production model may be overcome by studying ratios, e.g. + A)/(p + p)^ at the 
same center-of-mass energy since the dominance of gg processes means that the {p + A)/{p + p) ratio is, to a good 
approximation, the ratio of the gluon distribution in the nucleus relative to the gluon distribution in the proton. We 
have chosen to use the CEM because it allows predictions of the total cross section and the pT-integrated rapidity 
distributions where nuclear effects are more prominent. Measuring the J/ip and T ratios simultaneously also provides 
a means of determining the scale evolution of the nuclear gluon distribution at relatively large, perturbative scales if 
shadowing is the only cold nuclear matter effect in p + A and d+A collisions. 

At fixed-target energies, the xp dependence clearly shows that shadowing is not the only contribution to the J ftp 
nuclear dependence as a function of Xp [H, 111]. Indeed, the characteristic decrease of a{xF) for xp > 0.25 cannot 
be explained by shadowing alone [H. In fact, the data so far suggest approximate scaling with xf, not the target 
momentum fraction X2 [27'], indicating the possible importance of higher- twist effects \2§\ . The preliminary PHENIX 
data show an increasing suppression at forward rapidity ^29.] , similar to that seen in fixed-target experiments at large 

XF- 

Effects we have not considered here which may result in xf rather than X2 scaling and affect the high xf region are 
energy loss in cold matter and intrinsic charm, both discussed extensively in Ref. [Ij] . We do not consider these effects 
here because, at heavy-ion colliders, the relationship between xf, rapidity, and .^s„„ suggests that this interesting 
Xf region is pushed to far forward rapidities. The onset of initial-state energy loss should, in fact, appear at higher 
Xf 8ut larger ^s„„ if it depends on the momentum fraction xi. Figure [T] shows the relationship between xf and 
y in the center-of-mass frame for M = 4 GeV and 20 < ^s„„ < 14000 GeV. Since xf = {2mT / ^s„„ ) sinh y , the 
large center-of-mass energies at the LHC guarantees that the forward xf region will not be accessible in the central 
rapidity region of the LHC. Instead, the xf distribution becomes narrowly peaked with increasing energy while the 
rapidity distribution becomes broad and flat. At j/ = 5, the largest xf accessible (at the lowest ^s„„) is 0.081 for 
the J/ip and 0.25 for the T. The large xf region is therefore not probed by quarkonium production in \y\ < 5. Thus 
shadowing and absorption are likely the most important cold nuclear matter effects at the LHC. 

To implement nuclear absorption on quarkonium production in p + A and d-|-^ collisions, the production cross 
section is weighted by the survival probability, S(i'^, so that 

- s, z') = cxp |- dz"pA{b- s, z")ag,iz" - z')| (7) 

where z' is the longitudinal production point, as in Eq. ([3]), and z" is the point at which the state is absorbed. The 
nucleon absorption cross section, cr^g, typically depends on the spatial location at which the state is produced and 
how far it travels through the medium. 

We could also consider absorption by comover interactions but this cross section is typically smaller than the nuclear 
absorption cross section. In addition, in A -I- A collisions, the higher temperatures and larger particle densities would 
rule out hadronic comovers in the early stages. Therefore, we do not consider hadronic comovers as a cold matter 
effect in this paper. 

If absorption alone is active, i.e. Sp s{A, x, /i^, r, z) = 1, then an effective minimum bias A dependence is obtained 
after integrating Eqs. ([IJ and (O over the spatial coordinates. If S^}'^ — 1 also, (jpA ~ Aapp without any cold nuclear 
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matter effects. (Note that for gg-dominated processes, such as quarkonium production, the relationship would be 
exact. When qq or qq' interactions dominate, as in gauge boson production, the different relative proton and neutron 
numbers make the above relationship approximate.) If S'p g(A, x, fi^ ,r, z) = 1 and 7^ 1, apA = A^app where the 
exponent a can be related to the absorption cross section, as studied in detail for J/ip and tp' production by NA50 
[30|. For a constant cr^g with a sharp surface spherical nucleus of density pA — pq6{Ra -~ b), it can be shown that 

■rC 



■'0 



where rg = 1.2 fm [sT'l. The relationship between a and cr^g is less straightforward in more realistic geometries. 

The NA50 [S^l and E866 [25|| experiments measured a non-negligible difference in the effective J/-0 and ip' absorption 
cross sections at y^s^ = 23 — 29 GeV and ^/s^ = 38.8 GeV respectively. In addition, the difference between a'^^f 

and (T^i^g seems to decrease with ^/s„„. The NA50 collaboration measured Act = ct^j^j, — a'^l^ = 4.2 ± 1 mb at 400 
GeV and 2.8 ± 0.5 mb at 450 GeV At xp 0, the E866 resuhs imply Aa = aj/^ - a^, < 0.2 or, using Eq. ©, 
Act < 1.6 mb 2^]. This suggests that absorption is a final-state effect since an initial-state effect such as shadowing 
would not discriminate between the asymptotic J/tp and tjj' final states. Comparing the effective absorption cross 
sections determined at central rapidities from the CERN SPS to RHIC, absorption seems to decrease with energy 

M- 

Fewer T p + A data are available. The E772 experiment (3^ measured the A dependence of the three S states and 
found a reduced A dependence relative to J/ ip absorption. The A dependence of the three S states was indistinguish- 
able within the uncertainties. No T A dependence was presented by the E866 collaboration. The STAR d-|-Au/p-|-p 
ratio from RHIC suggests that, within large uncertainties, the T A dependence is linear [s^l and production is not 
significantly suppressed. Thus absorption seems to be weaker overall for T production but there is not clear indication 
so far of how much weaker it is or whether it has the same energy dependence as the J/ip. 

If conventional shadowing parameterizations, such as the ones used in this paper, are included, the effective J/ip 
absorption cross section may seem to decrease with energy due to the increased effect of shadowing at low x. A 
decrease in absorption concurrent with increased shadowing as ^s„„ increases seems to approximately hold, even 
without shadowing, at fixed-target energies (s^. Such a decrease is consistent with the J ftp traversing the nucleus as 
a color singlet. If the nuclear crossing time is shorter than the J/ip formation time, the effective absorption decreases 
with ^s„„ as an ever smaller state passes through the target. 

If the effective absorption cross section indeed decreases with energy, then absorption should be a relatively small 
contribution to the total A dependence at the LHC. This prediction is easy to check: if absorption is negligible, the 
J/ip and ip' {jp + A)/{p + p) ratios should depend only on shadowing and should thus be equivalent. The yield is 
then related to the ratio of the nuclear to proton gluon densities since gg fusion dominates quarkonium production at 
these energies. In this work, we have assumed that absorption is negligible so that the {p + A)/{p + p), {A+A)/{p + p) 
and {A + A)/{j) + p) J/^ and T ratios presented here are the same for all charmonium and bottomonium states 
respectively. 

If both the p + A and p+p data are taken at the same ^s„„ , the same x values of the gluon densities will be probed 
in the nucleus and in the proton. Such same energy comparison runs would be an excellent probe of the nuclear gluon 
distributions because 

°^ A cpi„ ..2\ ■ y^' 



jpp 



However, if the p -\- A and p + p data are recorded at different energies (and x values), the extraction of the nuclear 
gluon density is less straightforward since 

^ A CPl^ ..9.\ ■ \^^) 



Jpp 



(yi) A f^ix,p^) 



In both cases, the p^-integrated ratios provide an additional uncertainty because the scale evolution of the gluon 
density is not well known but is expected to be strong [sB-CSj. However, the quarkonium px distribution is steeply 
falling for pT > m so that the pT-integrated ratios ratios are a good representation of p,^ — {rriT)'^- 

The scale evolution of the gluon densities can be probed in part by relative studies of low pT or pT-integrated J/ip 
{m^ = 3.097 GeV) and T(IS') (mx(is) = 9.46 GeV) production. To more precisely obtain the scale evolution of 
shadowing, it would be preferable to bin the J/ip and T(IS') (p -I- A) /{p + p) ratios in px- One must be careful in the 
interpretation of such ratios, particularly at pT < m, since, at fixed-target energies, the pr-dependent (p + A) /{p + p) 



epender 

ratios show that the J/ip and T px distributions are broader in p+A than in p+p interactions [4lll4^ . This broadening 
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has been attributed to intrinsic parton px kicks accrued by the interacting parton as it traverses the nucleus before 
interacting [i^, . The magnitude of the average px kick increases with A so that the pT-dependent {p + A)/{p + p) 
ratio is less than unity at low pT and increases above one with increasing pT- This effect is important at low center- 
of-mass energies where the average pt of the produced quarkonium state is not large. By LHC energies, while the pt 
kick may be rather small relative to [p]^) , it may still affect the analysis of shadowing effects in pT-binned ratios but 
not in pT-integrated ratios. We focus on the pr-integrated results here and will present pT-dependent calculations 
elsewhere. 





CT''"^/nucleon pair (fJ,h) 




System y^s^ (TeV) 


J/'4^ Xcl Xc2 1p' 




P + P 14 


32.9 31.8 52.5 7.43 


3.15 0.055 


P + P 10 


26.8 26.0 43.3 6.06 


2.57 0.044 


P + P 9.9 
p+0 9.9 


26.6 25.8 42.6 6.02 
23.8 23.0 38.0 5.37 


2.55 0.044 
36.5 0.632 


P + P 9.39 
P+At 9.39 


25.8 25.0 41.3 5.83 
22.0 21.2 35.1 4.96 


2.48 0.043 
84.1 1.46 


P + P 9.27 
p+Kr 9.27 


25.6 24.8 40.9 5.79 
20.9 20.2 33.4 4.73 


2.46 0.043 
168.4 2.92 


P + P 9 
p+Sn 9 


25.2 24.4 40.2 5.69 
20.2 19.6 32.3 4.56 


2.41 0.042 
230.4 3.99 


P + P 8.8 
p+Pb 8.8 


25.0 24.2 39.9 5.65 
19.5 18.9 31.1 4.40 


2.40 0.042 
388.8 6.75 


P + P 7 
P+0 7 
d+O 7 
O+O 7 


21.8 21.1 34.9 4.93 
19.5 19.0 31.3 4.42 

19.5 19.0 31.3 4.42 

17.6 17.0 28.1 3.98 


2.09 0.036 
30.0 0.520 
60.0 1.04 
432.4 7.51 


P + P 6.64 
d+Ar 6.64 


21.2 20.5 33.8 4.78 
18.1 17.5 28.9 4.09 


2.02 0.035 
138.5 2.39 


P + P 6.48 
d+Kr 6.48 


20.9 20.2 33.3 4.71 
17.2 16.6 28.0 3.95 


2.00 0.034 
281.3 4.85 


P + P 6.41 
d+Sn 6.41 


20.7 20.1 33.1 4.68 

16.8 16.2 26.8 3.78 


1.98 0.034 
378.3 6.52 


P + P 6.3 
p+Ar 6.3 
Ar+Ar 6.3 


20.5 19.9 32.8 4.63 

17.6 17.0 28.1 3.97 
15.0 14.5 23.9 3.38 


1.97 0.034 
67.3 1.17 
2300 40.0 


P + P 6.22 
d+Pb 6.22 


20.4 19.7 32.5 4.60 
16.0 15.5 25.6 3.62 


1.95 0.34 
637.3 10.98 


P + P 6.14 
p+Kr 6.14 
Kr+Kr 6.14 


20.2 19.6 32.3 4.56 

16.6 16.1 26.6 3.76 

13.7 13.2 21.8 3.08 


1.94 0.034 
134.0 2.32 
9245 160.6 


P + P 5.84 
p+Sn 5.84 
Sn+Sn 5.84 


19.6 19.0 31.3 4.42 
15.9 15.4 25.4 3.59 
12.8 12.4 20.4 2.89 


1.88 0.033 
181.3 3.14 
17391 302.0 


P + P 5.5 
p+Pb 5.5 
Pb+Pb 5.5 


18.9 18.3 30.2 4.26 
14.9 14.4 23.8 3.37 
11.7 11.3 18.7 2.64 


1.81 0.032 
297.6 5.16 

48500 842 



TABLE II: The direct cross section per nucleon pair (central columns) and the dilepton yield per nucleon multiplied by AB. 
The results are given for the MRST PDFs [l3l with = 1.2 GeV, fiF ^ fJ-R = 2mT. 

As an example of the possible cross sections for quarkonium production at the LHC, we present the total cross 
sections in p+p, p + A, d+A and A + A collisions at the relevant energies. To illustrate the effects of shadowing on the 
total cross section, calculated to next-to-leading order in the CEM [Tgli we use the EKS98 parameterization |35l. [36|. 
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For each possible maximum p + A, d+A and A + A center-of-mass energy, we also give the p + p cross section at that 
same energy. In addition, for the A + A center-of-mass energies, we also give the p + p and p + A cross sections at 
that energy. The results are given in Tables |lT] and IIIII The central columns are the direct cross sections per nucleon 
pair for all charmonium and bottomonium states. The effects are largest for charmonium (lower x and than the T 
states) and for the heaviest nuclei (lowest energies - highest x - but largest A). On the right-hand side of the tables, 
the inclusive (direct plus feed down) cross sections are multiplied by the dilepton decay branching ratios. They are 
also multiplied by AB to obtain the minimum bias total cross sections. 

The approximate A dependence of the total cross section relative to the p + p cross section at the same center-of- 
mass energy, assuming no other cold matter effects, can be obtained from the (AB)" parameterization so that, per 
nucleon, 



a{p + A/p + p) 
a{A + B/p + p) 



In[//(4,A^^)//|(X2,M^)] 
\nA 

ln[//(x;,^^)4^(x^,A.^)/(/P(a:i,/i^)4n^2,A^^))] 
\n{AB) 



(11) 
(12) 



where = X2 and x'l ~ xi if the center-of-mass energies are the same for the two systems. For J/^p production in 
pPb and Pb-I-Pb collisions relative to p + p collisions at 5.5 TeV, a ~ 0.76 and 0.52 respectively. In the case of T 
production, we have a ^ 0.88 and 0.76 respectively. 

The p + p rapidity distributions for J/ip and T production at y/s^^^ — 5.5 and 14 TeV are compared to RHIC 
distributions at y/s^ = 200 and 500 GcV in Fig. [51 The LHC distributions are relatively constant over a range 
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-1 




-1 



3> 



-4 



FIG. 2: (Color online) The J/tp (left-hand side) and T (right-hand side) rapidity distributions at ~ 200 (dotted), 500 

(dot-dashed), 5500 (dashed) and 14000 (solid) GeV calculated as in Refs. [12, [il. The kinks in the J/ip distributions at LHC 
energies are the point where x < 10~^. Since the T factorization scale is larger, the T rapidity distributions are smoother. 
Note the different scales on the y-axes. 



of 5 or more units of rapidity, demonstrating that the cross sections are high enough to obtain good statistics for 
quarkonium states, even for forward production and detection, provided that the decay leptons are of sufficiently high 
Pt to reach the detectors^. 



^ This will be more difficult for CMS and ATLAS than for ALICE since the minimum muon pT for detection in the large p+p experiments 
is 3.5 GeV/c. 
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(j'^^'/nucleon pair (/ib) 


Ba'^'^AB (pb) 


System ,/s.,., fTeV) 


T T' T" -VhilP) \h(2P) 




P + P 14 


0.43 0.27 0.16 0.89 0.69 


0.020 0.0074 0.0036 


p + p 10 


33 21 12 70 54 


016 0059 0028 


p + p 9.9 
p+O 9.9 


0.32 0.20 0.12 0.66 0.51 
30 1 9 Oil 62 48 


0.015 0.0055 0.0026 
23 082 040 


p + p 9.39 
p+Ar 9.39 


0.30 0.19 0.12 0.63 0.49 
0.28 0.17 0.11 0.57 0.44 


0.014 0.0052 0.0025 
0.53 0.19 0.092 


p + p 9.27 

r)4-Kr 9 97 


0.30 0.19 0.11 0.62 0.48 
27 1 7 10 55 43 


0.014 0.0052 0.0025 
1 06 39 n 1 q 

-L-VJU U.Oi? U.-Lc/ 


P + P 9 
P_(_Sji 9 


0.29 0.18 0.11 0.61 0.47 
26 16 099 53 42 


0.014 0.0050 0.0024 
1 46 53 26 


p + p 8.8 

n+Ph 8 8 


0.29 0.18 0.11 0.60 0.47 
25 16 097 52 41 


0.014 0.0059 0.0024 
2 51 96 45 


P + P 7 
n+O 7 

d+0 7 
O+O 7 


0.23 0.15 0.090 0.48 0.38 
22 14 085 46 36 
0.22 0.14 0.085 0.46 0.36 
0.21 0.13 0.081 0.44 0.34 


0.011 0.0043 0.0019 
17 061 029 
0.34 0.12 0.058 
2.57 0.97 0.46 


71 -1- 71 fi4 

d+Ar 6.64 


22 14 085 46 36 
0.20 0.13 0.079 0.42 0.33 


01 1 0038 0019 

yj.yjA.±. \j yj.yjyji-rj 

0.78 0.27 0.13 


n 4- n fi 48 
d+Kr 6.48 


22 14 083 45 35 
0.20 0.12 0.076 0.41 0.32 


01 0037 001 8 
1.57 0.56 0.28 


p + p 6.41 
H+Sn 6 41 


0.21 0.14 0.082 0.44 0.35 
19 12 074 40 34 


0.010 0.0036 0.0018 
2 34 77 n 41 


p + p 6.3 
p+Ar 6.3 
Ar+Ar 6 3 


0.21 0.14 0.082 0.44 0.34 
0.20 0.12 0.075 0.41 0.32 
18 12 070 38 29 


0.010 0.0038 0.0018 
0.37 0.13 0.065 
13 8 5 29 2 43 

XtJ.O iJ.iJtJ i.^O 


p + p 6.22 
d+Pb 6.22 


0.21 0.13 0.080 0.43 0.34 
0.18 0.12 0.071 0.38 0.30 


0.010 0.0035 0.0017 
3.68 1.31 0.65 


n + n 6 14 
p+Kr 6.14 
Kr+Kr 6.14 


21 13 080 43 33 
0.19 0.12 0.072 0.39 0.30 
0.17 0.11 0.066 0.35 0.28 


0099 0038 001 7 
0.75 0.27 0.13 
57.4 21.8 10.1 


p + p 5.84 
p+Sn 5.84 
Sn+Sn 5.84 


0.20 0.12 0.076 0.41 0.32 
0.18 0.11 0.068 0.37 0.29 
0.16 0.10 0.062 0.33 0.26 


0.0094 0.0035 0.0017 
1.01 0.36 0.18 
108.1 41.3 19.0 


p + p 5.5 
p+Pb 5.5 
Pb+Pb 5.5 


0.19 0.12 0.070 0.39 0.30 
0.17 0.11 0.064 0.34 0.27 
0.15 0.094 0.057 0.31 0.24 


0.0090 0.0029 0.0016 
1.65 0.60 0.29 
304 116.1 53.5 



TABLE III: The direct cross section per nucleon pair (central columns) and the dilepton yield per nucleon multiphed by AB. 
The results are given for the MRST PDFs [3 with = 4.75 GeV, fip ^ fiR = rriT. 

II. COLD NUCLEAR MATTER EFFECTS 

In this section, we describe the cold nuclear matter effect of initial-state shadowing on J/iJj and T production at 
the LHC. We first discuss the shadowing parameterizations used in our calculations. We then show the effect of 
shadowing on the rapidity distributions in the p + A, d+A and A + A collisions available at the LHC. Finally, we 
discuss the collision centrality dependence on a simple model of inhomogeneous shadowing where the effect depends 
on the path length of the parton through the nucleus. 
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A. Shadowing parameterizations 



We use several parameterizations of the nuclear modifications in the parton densities to probe the possible range of 
gluon shadowing effects: EKS98 nDSg [13, HKN (sl, EPS08 ^ and EPS09 [40]. AU sets involve fits to data, 

typically nDIS data with additional constraints from other observables such as Drell-Yan dimuon production. Since 
these provide no direct constraint on the nuclear gluon density, it is obtained through fits to the dependence of 
the nuclear structure function, F^, and momentum conservation. The useful perturbative /Lt^ range of the nDIS data 
is rather limited since these data are only available at fixed-target energies. Thus the reach in momentum fraction, 
X, is also limited and there is little available data for x < 10~^ at perturbative values of /i^. This situation is likely 
not to improve until an eA collider is constructed p^. 

The EKS98 parameterization, by Eskola and collaborators, available for A > 2, is a leading order fit using the 
GRV LO proton parton densities as a baseline [H, [1^. The kinematic range is 2.25 < /^eks98 

< 10^ GeV2 

and 10~^ < X < 1. deFlorian and Sassot produced the nDS and nDSg parameterizations [s^l at both leading and 
next-to-leading order for 4 < A < 208. The weak gluon shadowing of the nDS parameterization appears to be 



ruled out by the rapidity dependence of J/ijj production at RHIC [47[. The stronger gluon shadowing of nDSg is 
used here. Calculations with the nDS parameterization predict negligible shadowing effects. The kinematic reach 
in X is the same as EKS98 while the /x^ range is larger, 1 < /i^Dgg < 10^ GeV^. Hirai and collaborators produced 
the leading order HKN parameterization by fitting parton densities for protons, deuterons and 16 heavier nuclei, 
typically those most commonly used in nDIS experiments. If a particular value of A needed for our calculations is not 
included, a set with a similar value of A is substituted. The HKN parameterization goes lower in x than the other 
parameterizations, 10~^ < a; < 1, and higher in scale, 1 < /^hkn < GeV^. The EPS08 parameterization, a fit 
by Eskola and collaborators that includes the BRAHMS d-|-Au data on forward rapidity hadron production at RHIC 
pSf . is designed to maximize the possible gluon shadowing^. The EPS08 x range is the same as EKS98, 10^^ < a; < 1, 
while the /i^ range was extended, 1.96 < /igpsos — GeV^. Very recently, the EPS09 [i^l parameterization, which 
excludes the BRAHMS data from the fits, was introduced. The EPS09 parameterization includes uncertainties on the 
global analyses, both at LO and NLO, by varying one of the 15 fit parameters within its extremes while holding the 
others fixed. The upper and lower bounds on EPS09 shadowing are obtained by adding the resulting uncertainties in 
quadrature The EPS09 central LO results are in quite good agreement with the older EKS98 parameterization 
while the maximum possible gluon shadowing effect resulting from their uncertainty analysis is similar to the EPS08 
gluon ratio. The minimal amount of gluon shadowing is nearly negligible, similar to nDS [37| and even leaves room 
for some antishadowing in light ions. We present the central EPS09 ratio as well as the ratios corresponding to the 
maximum and minimum range of the shadowing effect, obtained by adding the relative differences in quadrature, as 
prescribed in Ref. [40|. For computational convenience, we use the LO version of the nPDF parameterizations since 
the NLO GEM calculations give similar shadowing results (H^]. This is to be expected since, even though the LO and 
NLO values of the cross section and the shadowing parameterization are different, when convoluted, they give the 
same ratios by design, see e.g. Ref. (371] . 

While the x values probed at midrapidity are w 10^** for the J/ip and w 10^^ for the T, well within the x range 
of the parameterizations, this is not necessarily the case away from midrapidity. At the largest values of ^s„„, x 
values lower than the minimum valid x of the parameterization may be reached within the rapidity range of the LHG 
detectors. In these cases, the shadowing parameterizations are unconstrained by data. However, when x < 10^® 
the EKS98, nDSg, EPS08 and EPS09 parameterizations return the value of the shadowing ratio at x = 10~^. The 
minimum x value, 10~^, for the HKN parameterization is small enough that this minimum is not reached, even for 
the highest energies. 

The ratios of the nuclear gluon densities relative to the gluon density in the proton are shown in Fig. [3] for four 
different ion species available at the LHC: A — O, Ar, Sn and Pb. The calculations for A — Kr, an alternative 
intermediate mass ion species, are not shown. Results for scales appropriate for J/tp, Fig. [3Ja), and T, Fig. [H^b), 
production illustrate the scale dependence of the parameterizations. The scales correspond to those used in the 
calculations of the cross sections in Tables |ll] and IIIII with /i = 2mc for charm and nib for bottom respectively. If a 
lower scale, /i = rric, is used for charm, the shadowing effect is stronger since /i^ is then closer to the minimum scale of 
the parameterization. Note that in all cases shadowing increases with decreasing x and increasing A while decreasing 
with scale, /x, as seen by comparing Fig. (HJa) and (b). For example, the EKS98, nDSg and HKN ratios appear to be 
approximately independent of x for x < 10~^ at the J/i/j scale but not at the T scale. 

The EKS98, EPS08 and EPS09 parameterizations (solid and dotted curves and solid curves with symbols respec- 



^ It has been suggested that the BRAHMS data should not be used to calculate gluon shadowing effects and that the strong shadowing 
of EPS08 violates a unitarity bound at the minimum scale [49| . 
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FIG. 3: (Color online) The LO shadowing parameterizations for J/ip (a) and T (b) scales for O (upper left), Ar (upper right), 
Sn (lower left) and Pb (lower right) nuclei. The parameterizations are EKS98 (solid), nDSg (dashed), HKN (dot-dashed), 
EPS08 (dotted) and EPS09 (solid lines with symbols). Note that the lower limit on the y-axis is changed for Sn and Pb on the 
left-hand side. 



lively) exhibit large antishadowing, > 1, in the region 0.02 < a: < 0.2 — 0.3, becoming more pronounced for larger 
A. The nDSg parameterization (dashed curves) show very weak antishadowing around x ~ 0.1. At x < 10~^, the 
nDSg ratios are weakest for A = O and Ar, similar to HKN for A = Sn and compatible with EKS98 for x < 10~^. 
The HKN parameterization (dot-dashed curves), on the other hand, is similar to EKS98 for A = O but has a weak 
A dependence so that HKN shadowing is the weakest at low x and large A. The EPS08 parameterization is similar 
to EKS98 for x > 0.01 but exhibits stronger antishadowing at large A. It also has the strongest shadowing at low x 
since the \ow-pT forward-rapidity BRAHMS data was included in the fit. The scale dependence of nDSg and HKN 
appears to be weaker than EKS98. The EPS09 band is obtained by calculating the deviations from the central value 
for the 15 parameter variations on either side of the central set and adding them in quadrature. The range of the 
LO EPS09 uncertainty band encompasses all other shadowing ratios, similar to EPS08 for the maximum effect and 
even leading to antishadowing for lighter ions. (The central ratio is shown with circular symbols on the solid curve 
while the bounds include diamond symbols.) For smaller nuclei, the upper edge of the EPS09 uncertainty (minimal 
shadowing effect) gives a bound above unity for S^. 

All the parameterizations increase at large x with 5^ > 1 for a; > 0.1 (HKN and nDSg) and x > 0.7 (EKS98 and 
EPS08). The rise in the HKN parameterization is steepest and occurs at the lowest x, beginning at the x value of the 
antishadowing peak in the EKS98 and EPS08 ratios. This high x region will not be explored by the LHC detectors 
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since it is only reached at rapidities outside their acceptance. 

Finally, we note that since our p + A calculations assume the ion beam travels in the negative z direction, low x 
corresponds to large forward rapidity while high x corresponds to large backward rapidity. 
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FIG. 4: (Color online) The (p + A)/(p + p) ratios with both p + A and p + p collisions at the p + A energy in the equal-speed 
frame. No rapidity shift has been taken into account. The effect of shadowing on J/tp (a, upper 4 panels) and T (b, lower 
4 panels) production is shown. Each panel displays the production ratios for p+O at .^s„„ = 9.9 TeV (upper left), p+Ar 
at ./s^ = 9.39 TeV (upper right), p+Sn at ^/s^ = 9 TeV (lower left) and p+Pb at = 8.8 TeV (lower right). The 

calculations are with CTEQ6 and employ the EKS98 (solid), nDSg (dashed), HKN (dot-dashed) and EPS09 (solid curves with 
smbols) shadowing parameterizations. 



B. Rapidity dependence 

We now show predictions of the J/ip and T production ratios as a function of rapidity for cold nuclear matter, 
CNM, effects at the LHC. li h + A data (where h = p ot d) can be taken at the same energy as the p + p and/or A + A 
data, as at RHIC, it is easier to make comparisons. However, the setup of the LHC makes this ideal situation more 
difficult. At the nominal injection energy, the proton beam has an energy of 7 TeV while the nuclear beam energy 
per nucleon is lower by the nuclear charge-to-mass ratio, Z/A. To make a p + p comparison, if we are not to rely 
on calculations extrapolated to lower energy, the p + p collisions have to be run at the p + A oi A + A per nucleon 
energies. For the proton and ion beam energies to be the same, the proton beam must then circulate at lower than 
optimal energy, decreasing the luminosity. Since sustained low energy p + p runs are unlikely in early LHC running, 
especially for sufficiently accurate quarkonium data as a function of rapidity, it may be necessary to rely on higher 
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energy p + p reference data'^. However, there is a catch. In p + A colhsions where a 7 TeV proton beam colhdes with 
a 7{Z/A) TeV per nucleon ion beam, the so-called equal-speed or equal-rigidity frame, the center-of-mass rapidity is 
not fixed at y = but displaced by Ay^^. In pPh collisions, the shift can be nearly 0.5 units, an important difference, 
see Table |T] for the magnitude of the possible shifts. To minimize the rapidity shift and to bring the hA comparison 
energy closer to that of the A + A energy, d+A collisions may be desirable since = 3.5 TeV per nucleon relative to 
Epb = 2.75 TeV, see Table HI Since d+A collisions require a second ion source, this may not be realized in the short 
term. 

We thus study several different possibilities for determining cold nuclear matter effects on nucleus-nucleus collisions 
at the LHC. We go from ideal to more realistic scenarios. We first show the {p + A)/{p + p) per nucleon ratios at the 
same per nucleon center-of-mass energy for both systems, assuming the appropriate p+p energies are available, Figs.|4] 
and[SJ In the case where p-\- A and p + p interactions are compared at the A + A energy, we assume zero rapidity gap, 
Aj/cm — 0) between the colliding beams since the proton beam energy is reduced to match that of the nucleus. In the 
more likely scenario. Figs. [5] and [71 the p + A data will be taken in the equal-speed frame at a higher energy than 
the A + A collisions. Therefore, we next show the {p + A)/{p + p) per nucleon ratios with respect to p + p collisions 
at yfs = 14 TeV with Ay^^ — for p + A collisions both in the equal-speed frame and at the A + A center-of-mass 
energy. The final p + A calculations shown are the most realistic: the p + A cross section in the equal-speed frame with 
finite AyP^ is shown relative to the p + p cross section at 14 TeV. In this case, Fig.[5J the numerator and denominator 
are calculated with different energies and different center-of-mass rapidities. Next, the {d+A)/{p + p) per nucleon 
ratios are presented for two cases: with the d+A and p + p collisions at the d+A center-of-mass energy and with d+A 
collisions in the equal-speed frame with Ay'^^ ^ and p + p collisions at 14 TeV. Finally, we present the baseline 
(A + A)/[p + p) per nucleon ratios with the p + p center-of-mass energy tuned to the A + A energy and at the nominal 
14 TeV p + p energy. In the case of symmetric p + p and A + A collisions, there is no rapidity gap. 

We present the rapidity dependence of p + A, d+A, and A -I- A collisions for A — 0, Ar, Sn and Pb relative to p + p 
collisions, both at the same energy as the nuclear system and at 14 TeV. The A + B/p + p ratios are shown for the 
EKS98, nDSg and HKN parameterizations. The EPS09 central ratios and the associated uncertainty bands are also 
shown. Since the EPS08 ratios are similar to the lower limit (strongest shadowing) of the EPS09 uncertainty band 
at small x, we do not show any further calculations with EPS08. We use the CTEQ6 parton densities to calculate 
the ratios shown in the next two sections. We have checked that the ratios with the MRST densities are essentially 
identical since the same mass and scale parameters are used in the two calculations*. Both J /tjj and T results are 
shown in each figure. To guide the reader and clarify the discussion, we first present a table of the figures with the 
center-of-mass energy of the A + B and p + p collisions and the rapidity shift. 
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number 



Cross Section Ratio 

{B + A)/{p + p) 



pp 

Sri 



O'BA 



(TeV) ^iiX; (TeV) ySjjl (TeV) y/I^ (TeV) 



Ay, 



BA 
cm 



B 



m 
m 
m 
m 



<^va{\/saa, y)/[Aapp{^SAA, y)\ 
<JpA{^/s^, y)/[Aapp{^/s^,y)] 
e^pAiy/s^, y)/[Aapp{y/s^, y)] 
apA{y/s^, {y - AyP^)))/[Aapp{y/8^, y)] 



^ SpA 


9.9 


9.39 


9 


8.8 





yfSAA 


7 


6.3 


5.84 


5.5 





14 TeV 


9.9 


9.39 
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8.8 





14 TeV 
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6.3 


5.84 


5.5 





14 TeV 


9.9 


9.39 


9 


8.8 


Ay?^ 



B ■ 



m 



(y&A{yfsAA, y)l\2Aapp{^SAA, y)\ 
Ay,^^))/[2A<7j,p(y5^ 



6.64 
6.64 



6.41 
6.41 



6.62 
6.62 



B = A 



m 

[12] 



GAA{s/s^,y)l\A Opp{yfs^,y)\ 

oaa{^saa, y)l{J^app{yfs^, y)] 



6.3 
6.3 



5.84 
5.84 



5.5 
5.5 



^ Although the startup LHC %> + v runs are at lower energies, it is not clear how much quarkonium data will be extracted during these 

runs. Thus we show our results relative to the maximum p + p energy of 14 TeV. 
^ The ratios would only differ if these parameters were changed. However, we leave them fixed since they were optimized to other data 

for a given set of parton densities. 
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TABLE IV: Summary of the contents of Figs. |4]-[l2]in Section [II B I Here B is the identity of the colUsion partner, B — p for 
p + A, d for and A for ^ + ^4 collisions. The value of the center-of-mass energy for p + p collisions used in the calculation 
of the baseline p + p cross section is given in the third column: ^Spp — y^SpA for p + A; ^/sdA for d+A; ^/saa for A + A; and 
14 TeV for maximum energy p + p collisions. The center-of-mass energy for the B + A cross sections are given in the next four 
columns. Finally, whether or not the rapidity shift is included is indicated in the last column. The value of Ay, 
Table U] Note that all ratios are given for the per nucleon B + A cross section. 



^ IS given m 



The {p + A)/{p + p) ratios with equal p + A and p + p center-of-mass energies, shown in Figs. S] and [SJ illustrate the 
direct shadowing effect. The ratios are given both at the energy in the equal-speed frame, the likely ^s„„ for p + A 
collisions (Fig. S]), and at the same ^s„„ as the corresponding A + A collisions (Fig. [5]). The results are shown for all 
shadowing parameterizations. The nuclear beam is assumed to be moving from positive to negative rapidity so that 
the smallest values of x probed in the nucleus are at large, positive y. 

The LHC could be run as either ap + A or a.n A+p collider. Since the ATLAS and CMS detectors are symmetric 
around y = with central muon detectors in the range \y\ < 2.4, ALICE is the only experiment that could benefit 
from running in both modes because their dimuon spectrometer covers —4 < y < —2.4 in these coordinates fsH . 
However, since ALICE has muon coverage in the largest y region, running in both modes could be an advantage for 
reconstructing the nuclear effects in quarkonium measurements, especially since the y distributions are rather flat 
over a broad rapidity range. The large rapidity rates are thus non- negligible. 
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FIG. 5: (Color online) The (p + A)/{p + p) ratios with both p + A and p + p collisions at the A + A center-of-mass energy. The 
effect of shadowing on J /ijj (a, upper 4 panels) and T (b, lower 4 panels) production is shown. Each set of panels displays the 
production ratios for p-j-O at ^s„„ = 7 TeV (upper left), p-j-Ar at ys„„ — 6.3 TeV (upper right), p-j-Sn at — 6.14 TeV 

(lower left) and p-j-Pb at .ys„„ = 5.5 TeV (lower right). The calculations are with CTEQ6 and employ the EKS98 (solid), 
uDSg (dashed), HKN (dot-dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. 
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The J /ip ratios are shown in the upper half of the figures while the T results are on the lower half. Since shadowing is 
an initial-state effect, the same ratios would also be expected for the Xc and V'' on the left and the higher T states (T', 
T", Xf,(lP) and Xbi'^P)) on the right. The ratios in Figs. |3] and [5] are stretched mirror images of the gluon shadowing 
ratios in Fig. [31 The lowest x values are probed by the lightest nuclei since the center of mass energy is higher for 
nuclei with Z/A ~ 0.5 than heavier, neutron-rich nuclei with lower Z/A. The differences in the shadowing ratios for 
a given parameterization are greatest at large negative y where x is largest. As A increases and y/s^^ decreases, the 
antishadowing peak moves closer to midrapidity (less negative y). Increasing the scale from that appropriate for the 
J /ip to that for the T also moves the antishadowing peak closer to y — Q. For example, the EKS98 antishadowing 
peak is fully visible for T production, ocurring at y '-^ —3.5, while they only appear at y < —5 for the J/ip. As 
/i^ increases, the differences in the EPS09 sets becomes more pronounced at large x, leading to the more irregular 
shapes of the upper and lower limits of the EPS09 uncertainty range at negative rapidity. Note that the central ratio 
is smooth. Thus, the results in Figs. |4] and [5] suggest that by running the LHC in both p + A and A + p modes the 
modification of the nuclear gluon parton density could be traced out over a wide x range, taking advantage of the 
ALICE muon coverage. 

Finally, we note that at y = 6, corresponding to a; < 10~^, the EKS98 and nDSg shadowing ratios are outside their 
range of validity. This is also near the region where DGLAP evolution of the parton densities is likely to break down. 
Nonlinear evolution of the proton parton densities is expected at sufficiently small x. The onset of these nonlinearities 
is predicted to be at larger x for nuclei. However, it is not obvious that nonlinear parton evolution automatically 
leads to a reduction of the small x gluon density even though the nonlinear term in the gluon evolution has a negative 
sign (H^]. The behavior of the low x gluon density cannot be determined without a complete re-evaluation of all the 
parton densities since the sea quark evolution is coupled to that of the gluon and overall momentum conservation 
must be maintained along with the integrity of the global fit. See Ref. 53] for details of modified parton densities 
based on nonlinear DGLAP evolution and Refs. [5J-l56[ for a discussion of the possible effect on charm production at 
the LHC. 

Since it is more likely that the best p + p reference data will be at y^s — 14 TeV or 10 TeV for the initial LHC run. 
Figs. [6] and [7] show the {p + A)/{p + p) ratios with the p + p reference at 14 TeV. The magnitude of the two ratios 
(for p + A collisions in the equal-speed frame. Fig. [SJ and at the same energy as the corresponding A + A collisions, 
Fig. [7]) is due to the difference in ^/s^ relative to 14 TeV. The p + A ratios with the p + A center-of-mass energies 
equal to those oi A + A collisions are lower. The rapidity distributions narrow while their magnitudes are reduced 
with decreasing ^s„„. Thus the {p + A)/{p + p) ratios without shadowing decrease steadily from 9.9 to 5.5 TeV while 
the narrowing of the ratios becomes more pronounced. 

The symmetric solid curves in Figs. [S] and [7] are the {p + A)/{p + p) ratios without shadowing. Shadowing results in 
asymmetric ratios but since the p + A phase space is narrower than that of 14 TeV p + p collisions, the ratios in these 
figures turn over and drop to zero at large \y\. The narrower phase space has a bigger effect on the T production 
ratios since the full T rapidity range is within \y\ < 6 while the J/^ y distribution is broader. The antishadowing 
peak is lowered and broadened when dividing by the 14 TeV p + p rapidity distribution and is only really apparent for 
the EKS98 and EPS09 parameterizations. The maximum shadowing allowed by EPS09 shows the most asymmetric 
curvature, especially for J/tp- The EPS09 ratios suggest that the effect could either be large, as suggested by the 
EPS08 analysis, or small enough to be effectively indistinguishable from no shadowing. It will thus be harder to 
differentiate between shadowing parameterizations when employing the higher energy p + p reference. 

As discussed previously, there is an additional complication due to the rapidity shift of the p + A center of rapidity 
in the equal-speed frame. The shift increases with A as Z/A decreases, reducing the energy of the ion beam relative 
to the proton beam. This results in nearly half a unit rapidity shift in pPb collisions, as shown in the center part 
of Table [H labeled p + A. The {p + A)/{p + p) ratios including the rapidity shift and the maximum energy p + p 
reference are shown in Fig. [5] Note that only the p + A results in the equal-speed frame are shown. Since the proton 
beam momentum in p -I- ^ collisions at the A + A center-of-mass energy must be the same as that of the ion beam, 
^Ucm — 0- The p + A rapidity distribution is given a positive shift, to the right, since the proton beam, at higher y, 
is assumed to come from the left and move to the right. Thus, at large negative y, the ratios are lower than in Figs. [6] 
and [7] and are flatter as a function of rapidity. While the nuclear effects on the parton densities are most difficult to 
disentangle here, this scenario is the most realistic. (It may be possible to eliminate or reduce the effect of the rapidity 
shift by employing different rapidity cuts to compare distributions.) If the LHC is run with the proton and ion beam 
directions reversed, the antishadowing peak may be enhanced and the large positive rapidity ratios decreased. 

The effect of the rapidity shift is reduced if d-|-A collisions are run instead oi p + A collisions. The d-f ^ center- 
of-mass energy is closer to that oi A + A collisions since Z/A < 1 for the deuteron rather than equal to 1 as for 
protons. The ratios with the d+A and p + p collisions at the same center-of-mass energy per nucleon are shown in 
Fig. [9] (similar to Fig. [5] for p + A). They are like those in Fig. [5] with equal p + A and A + A center-of-mass energies 
since ^s„„ is similar for d+A and A + A collisions. Shadowing effects on the deuteron are assumed to be negligible. 

The results with a 14 TeV p + p reference and the small rapidity shift taken in account, see Table U for Ay'^^, 
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FIG. 6: The (p + A)/ (p+p) ratios with the p+p rapidity distributions calculated at ^/s = 14 TeV. While the p + A distributions 
are calculated in the equal-speed frame, no rapidity shift has been taken into account. The effect of shadowing on J/^ (a, 
upper 4 panels) and T (b, lower 4 panels) production is shown. Each set of panels displays the production ratios for p+O at 
= 9.9 TeV (upper left), p+Ar at ^/s^^ = 9.39 TeV (upper right), p+Sn at ^/sj^ = 9 TeV (lower left) and p+Pb at 
^/s'^ = 8.8 TeV (lower right), all calculated in the equal-speed frame. The calculations are with CTEQ6 and employ the 
EKS98 (solid), uDSg (dashed), HKN (dot-dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. The 
solid curve symmetric around y = is the {p + A)/{p + p) ratio without shadowing. 



are shown in Fig. [101 Recall that there is no rapidity shift for d-|-0 collisions since Zd/Ad = Zq/Aq = 0.5. Thus 
the equal-speed and center-of-rapidity frames coincide. In d-|-Pb collisions, since Ay;^^ < 0.06, the shift is negligible. 
Thus the d+A rapidity distributions relative to the 14 TeV p + p reference with the rapidity shift, shown in Fig. [TOl 
are similar to those in Fig. [7] with Ay^^^ = and the same ^s„„ m p + A and A + A collisions. Note, however, 
that the ratios in Fig. [TU] are somewhat closer to unity since the d+A center-of-mass energy is larger. Thus the more 
realistic d+A scenario shown in Fig. [TU] would be preferable for determining nuclear effects on the parton densities 
both because of the relatively similar center-of-mass energies and the smaller rapidity shift. 

We now extrapolate to A + A interactions to show the projected CNM effects from shadowing alone. The results 
for A + A collisions are shown in Figs. [TT] and [T^l The {A + A)/{p + p) ratio with both systems calculated at the 
A + A center-of-mass energy are shown in Fig. [TT] while the 14 TeV p + p reference is employed to obtain the ratios 
in Fig. [T^] The results in Fig. [TT] are essentially the convolutions of the {p + A)/{p + p) ratios (with the same ^/s^ 
for both systems and no rapidity shift) shown in Fig. [5] with their mirror image Ap/pp ratios. While the A + A ratios 
exhibit antishadowing peaks at y ~ ±(4 — 5), the {A + A)/{p + p) ratios are less than unity everywhere because the 
product of the {p + A)/(p + p) ratios at postive y and the {A + p)/ {p + p) ratios at negative y is always smaller than 
one, e.g. [{p + A)/{p + p)]y~5 ^ 0.6 — 0.75 while [{A + p)/{p + p)]yr^-5 ^ 1.2. Thus, when all ratios are calculated 
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FIG. 7: The {p + A)/{p + p) ratios with the p + p rapidity distributions calculated at y/s — 14 TeV. The p + p distributions 
are calculated at the A + A center-of-mass energy. The effect of shadowing on J/^p (a, upper 4 panels) and T (b, lower 4 
panels) production is shown. Each set of panels displays the production ratios for p+O at y^s^Tiv ~ ^ TeV (upper left), p+Ar 
at .ysjyjy = 6.3 TeV (upper right), p+Sn at .ys„„ = 6.14 TeV (lower left) and p+Pb at ^^s^„ ~ 5.5 TeV (lower right). The 
ratios are calculated at the A + A energy. The calculations are with CTEQ6 and employ the EKS98 (solid), uDSg (dashed), 
HKN (dot-dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. The solid curve symmetric around 
y = is the (p + j4)/(p + p) ratio without shadowing. 



at the A + A center-of-mass energy, assuming factorization of A + A collisions into a convolution of p + A and A+p 
collisions, 

A + A P + A A + p 

; \y~±5 = ; y~5 X ; \y 5 < 1 • (13) 

p + p p + p" p + p^ 

Calculations of color singlet J/-0 interactions in matter using the dipole approximation of the J/T/i-hadron cross 
section suggest that factorization is inapplicable due to the coherence of the interaction [s^l • These gluon saturation 
models assume the dominance of higher- twist effects enhanced by powers of A^/^ \ap + A interactions. If these models 
were valid, enhanced suppression of the J /ij) should set in at large rapidity. This does indeed seem to be the case at 
RHIC where 1.2 <y < 2.2 corresponds to 0.0045 > X2 > 0.0017 [i3|- However, the forward xp data at = 38 GeV 
[0.2 < xfO.S and 0.027 > X2 > 0.008] and 17 GeV [0.1 < xp < 0.35 and 0.13 > X2 > 0.07], in particular, exhibit 
the same trends as at RHIC [H^. These fixed-target X2 ranges lie in the transition region from antishadowing to 
shadowing (38 GeV) and in the antishadowing region (17 GeV), see Fig.[3l seemingly too large to be in the saturation 
region, especially at ^/s — 11 GeV. 

As is the case for the RHIC A + A calculations at ^s„„ — 200 GeV [11], there is typically more suppression 
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FIG. 8: (Color online) The {p + A)/(p + p) ratios with the p + p rapidity distributions calculated at ^/s = 14 TeV. The p + A 
rapidity distributions are calculated in the equal-speed frame with the rapidity shift taken into account. The effect of shadowing 
on J ftp (a, upper 4 panels) and T (b, lower 4 panels) production is shown. Each set of panels displays the production ratios 
for p+0 at y^s^Tiv" ~ 9-9 TeV (upper left), p+Ar at v/sTriv ~ 9-39 TeV (upper right), p+Sn at y^sTriv ~ 9 TeV (lower left) and 



p+Pb at 



= 8.8 TeV (lower right). The calculations are with CTEQ6 and employ the EKS98 (solid), nDSg (dashed). 



HKN (dot-dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. The upper solid curve at y > is the 
shifted {p + A)/{p + p) ratio without shadowing. 



predicted &i y — Q than at more forward and backward rapidities for all the shadowing parameterizations as well 
as for both J/-0 and T production. At RHIC, the A + A data are more suppressed at forward rapidity than at 
central rapidity, both in the minimum bias data as a function of rapidity and as a function of collision centrality, as 
quantified by the number of participant nucleons. Standard models of shadowing alone or shadowing with absorption 
by nucleons in cold nuclear matter or shadowing combined with dissociation in a quark-gluon plasma leads to str ong 
suppression at central rapidities. However, J/'0 regeneration by coalescence of c and c quarks in the medium l60l | 
is biased toward central rapidities and could lead to more suppression at forward rapidity relative to central rapidity 
since the rapidity distribution of J/ip production by coalescence is expected to be narrower than the initial J/ij] 
rapidity distribution [6^. Thus, with coalescence, there should be more suppression at forward y than at midrapidity. 
The same trend should hold at the LHC. Coalescence production of the J/-0 should be even more important than at 
RHIC since more cc pairs are created in a central Pb-|-Pb collision at y/s^ = 5.5 TeV. We can also expect that T 
production by coalescence may be similar to that expected for the J/tp at RHIC since the bb production cross section 
at the LHC will be similar to the cc production cross section at RHIC [13] . 

The (A + A)/{p + p) ratios with the 14 TeV p + p reference, shown in Fig. [121 a-i'e relatively flat. The dip around 
midrapidity has been washed out, except for the J/ip ratios calculated with the EKS98 and EPS09 (central and 
maximum shadowing) parameterizations where some indication remains. For comparison, the {A + A)/(p + p) ratios 
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FIG. 9: (Color online) The {d+A)/{p + p) ratios with both d+A and p + p collisions at the d+A energy in the equal-speed 
frame. No rapidity shift has been taken into account. The effect of shadowing on J/tp (a, upper 4 panels) and T (b, lower 4 
panels) production is shown. Each set of panels displays the production ratios for d+O at ys„„ = 7 TeV (upper left), d+Ar 
at = 6.64 TeV (upper right), d+Sn at ^/s^ = 6.41 TeV (lower left) and d+Pb at ^/s^ = 6.2 TeV (lower right). The 

calculations are with CTEQ6 and employ the EKS98 (solid), uDSg (dashed), HKN (dot-dashed) and EPS09 (solid curves with 
symbols) shadowing parameterizations. 



without shadowing are shown in the upper solid curves. 



C. Impact parameter dependence 



We now discuss the impact parameter dependence of quarkonium production at the LHC. Unfortunately, there is 
little relevant data on the spatial dependence of shadowing. Fermilab experiment E745 studied the spatial distribution 
of nuclear structure functions with vN interactions in emulsion. The presence of one or more dark tracks from slow 
protons is used to infer a more central interaction [6ll|. For events with no dark tracks, no shadowing is observed while, 
for events with dark tracks, shadowing is enhanced over spatially-independent measurements from other experiments. 
Unfortunately, this data is too limited to be used in a fit of the spatial dependence. 

The minimum bias shadowing we have discussed up to now is homogeneous, impact parameter-integrated shadowing. 
The impact parameter-dependent results shown in this section portray inhomogeneous shadowing. In central collisions, 
with small impact parameter, b, we can expect inhomogeneous shadowing to be stronger than the homogeneous result. 
In peripheral (large impact parameter) collisions, inhomogeneous effects are weaker than the homogeneous results but 
some shadowing is still present due to the overlapping tails of the density distributions. The stronger the homogeneous 
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FIG. 10: (Color online) The {d+A)/(p + p) ratios with the p + p distributions calculated at -y/s — 14 TeV and the d+A 
rapidity shift taken into account. The effect of shadowing on J/ip (a, upper 4 panels) and T (b, lower 4 panels) production is 
shown. Each set of panels dispalys the production ratios for d+O at ,ys^^ = 7 TeV (upper left), d+Ar at y^s^^iv" ~ 6.64 TeV 
(upper right), d+Sn at .^s„„ = 6.41 TeV (lower left) and d+Pb at ^/sj^ = 6.2 TeV (lower right). The calculations are with 
CTEQ6 and employ the EKS98 (soUd), nDSg (dashed), HKN (dot-dashed) and EPS09 (solid curves with symbols) shadowing 
parameterizations. The symmetric solid curve is the result without shadowing. 



shadowing, the larger the difference between the central and peripheral results. 

We assume that the shadowing is proportional to the parton path through the nucleus [s^l , 

S^pJA,x,Q\r,z) = l + N,{Si,{A,x,Q')-l)j^^^^ , (14) 

where Np is chosen to satisfy the normalization condition in Eq. ([S]). The integral over z in Eq. (|14p includes the 
material traversed by the incident nucleon. At large distances, s S> Ra, the nucleons behave as free particles, while 
in the center of the nucleus, the modifications are larger than the average value Sp. 

We calculate the nuclear suppression factor, Rab, for p + A, d+A and A + A collisions. The suppression factor is 
defined as the ratio [63| 

where daAB /dy and dopp/dy are the quarkonium rapidity distributions va. A + B and p + p collisions and Tab is the 
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FIG. 11: (Color online) The {A + A)/{p+p) ratios with both A + A and p + p collisions calculated at the A + A center-of-mass 
energy. The effect of shadowing on J/tp (a, upper 4 panels) and T (b, lower 4 panels) production is shown. Each set of panels 
displays the production ratios for 0+0 at ^s„„ = 7 TeV (upper left), Ar+Ar at .ys„„ = 6.3 TeV (upper right), Sn+Sn at 
.ys„„ = 6.14 TeV (lower left) and Pb+Pb at ,ys^^ — 5.5 TeV (lower right). The calculations are with CTEQ6 and employ 
the EKS98 (solid), nDSg (dashed), HKN (dot-dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. 



nuclear overlap function, 

TAsib) = J(fsdzdz'pA{s,z)pBi\b~s\,z'). (16) 

In p + A collisions, we assume that the proton has a negligible size, pa{s, z) — S{s)6{z) so that 2ab(6) collapses to 
the nuclear profile function Tsib) = J dz' pB{b, z'). The deuteron cannot be treated like a point particle since it is 
large and diffuse. We use the Hiilthen wave function [gjj to calculate the deuteron density distribution. However, we 
do not include shadowing effects on the deuteron. 

We show the p+Ph and d+Pb suppression factors as a function of impact parameter with ^/s^ — 8.8 TeV and 6.2 
TeV in the numerator and denominator in Figs. [T3l and [T4l respectively. We concentrate on the largest A ion, Pb, to 
maximize the relevant impact parameter range. The results in Fig. [13] are given for three values of rapidity: y = —4 
(backward rapidity, in the antishadowing region for T), y = (midrapidity) and y = 4 (forward rapidity, where fairly 
strong shadowing is expected). We present J/ip ratios on top and T ratios on the bottom. For comparison, the 
horizontal lines, centered around the average path length through the lead nucleus, b ^ (3/4)i?pb, show the impact 
parameter- integrated ratios in Fig. 21 The b dependence is strong, resulting in Rpph ^ 1 for b > Rph- Shadowing is 
stronger in central colisions than the average integrated value, as expected. Because the average decreases at forward 
rapidities while the spatial dependence is relatively unchanged, the strongest b dependence is seen for the most forward 
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FIG. 12: (Color online) The (A + A)/(p + p) ratios with the p + p rapidity distributions calculated at ^/s = 14 TeV. The 
effect of shadowing on J/?/) (a, upper 4 panels) and T (b, lower 4 panels) production is shown. Each set of panels displays the 
production ratios for O+O at .ys„„ = 7 TeV (upper left), Ar+Ar at ^/s^ = 6.3 TeV (upper right), Sn+Sn at ^/s^ ~ 6.14 
TeV (lower left) and Pb+Pb at ^/i]^ = 5.5 TeV (lower right). The calculations are with CTEQ6 and employ the EKS98 
(solid), nDSg (dashed), HKN (dot-dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. The upper 
solid curve is the {A + A)/{p + p) ratio without shadowing. 



rapidity value, y = 4. The nDSg, EKS98, and EPS09 shadowing ratios in Fig. |3] are very similar for lead nuclei. Thus 
their suppression ratios are also similar. Since only one nuclear density profile is involved in the calculation of RpA, 
the impact parameter dependence reflects that of Eq. [14] rather directly. 

A weaker impact parameter dependence is seen for d+Pb collisions in Fig. [T31 The overall shadowing effect is 
reduced since the energy is lower, ^/s^ — 6.2 TeV relative to 8.8 TeV for pPh collisions. In addition, shadowing 
persists to large values of impact parameter. In a heavy nucleus, the density is large and approximately constant except 
close to the surface, as expressed by the Woods-Saxon density distributions [l^. However, the diffuse wavefunction of 
the deuteron has a finite amplitude at surprisingly large distances. These lon g ta ils produce some remnant shadowing 
effect even at very large b, as seen in Fig. [TH See also the discussion in Ref. [63 ■ 

The d-|-Au results at RHIC have been presented as a function of the number of binary nucleon-nucleon collisions, 
-^coii(sjvjv ; ^) = Cinei(s„jv )Tab(&), rather than impact parameter itself, see Ref. |47| for details. The number of 
collisions is greatest for the most central collisions, 6 « 0, and decreases with increasing b. Since the inelastic nucleon- 
nucleon cross section, ainciis j^) , is energy dependent, the number of collisions increases with energy even though 
TAsib) does not. Thus Ncoi\{s^„;b) is significantly larger at the LHC than at RHIC for the same A + B system 
because of the considerable increase in (Tinoi(s„„) (from 42 mb at y/J^ = 200 GeV to 75 mb at = 5.5 TeV. 

The results for nucleus- nucleus collisions are presented as a function of the number of participant nucleons, iVpart, 
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FIG. 13: (Color online) The suppression factor Rpph at y — —4 (left), (center) and 4 (right) as a function of 6. The result is 
shown for J/tp (top) and T (bottom) in p+Pb relative to p + p colHsions at the same energy, y^s^^jv" ~ 8-^ TeV, and employ 
the EKS98 (solid), nDSg (dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. The horizontal lines 
show the impact-parameter integrated results. 



which depends on b as 

iVpart(fo) = / d^s [r^(s)(l-exp[-ai„ei(s„„)TB(|6-.^)]) 



+TB(|6-.s1)(l-exp[-ai„ei(s„„)r^(s)])J . (17) 

Large vahies of iVpart are obtained for small impact parameters with A^partl^ = 0) = 2 A for spherical nuclei. Small 
values of A^part occur in very peripheral collisions. Figure [TSl shows RAA{Npart) , at y ~ for the four A + A systems 
where the p + p and A + A rapidity distributions are calculated at the same center-of-mass energy. A similar pattern 
is observed for other values of y since the {A + A)/(jj + p) ratios are approximately independent of rapidity over a 
rather broad range. The {A + A)/{p + p) ratio at y = from Fig. [TT]is indicated by a horizontal line. Note that 
RAA{Npart) in Fig. [15] is equal to {A + A)/{p + p) in Fig. [TT] for Afpart(^ « Ra)- In smaU systems, i?AA (A'^part) is 
almost linear with more curvature appearing for larger collision systems. 

Since the p+p reference is not likely to be immediately available at the A + yl center-of-mass energy for Raa studies, 
Eq. (|15p . it may be preferable to study ratios of two quantities measured at the same energy in A + B collisions where 
B = p, d, or A. In this case, we utilize Rcp, the ratio of A + B cross sections in central relative to peripheral collisions, 

jl^^pl^yj = Tab [bp) daAB{hc)ldy ^^^^ 
TAB{bc) d(7AB{bp)/dy 

where be and bp correspond to the central and peripheral values of the impact parameter. Indeed, shadowing may best 
be probed by Rcp measurements in asymmetric systems since the most peripheral collisions are a good approximation 
to nucleon-nucleon collisions. The same rapidity shift is common to both central and peripheral collisions. We note, 
however, that an experimental measurement will not be able to define a precise impact parameter but will instead 
define impact parameter bins of finite width. Thus any comparison of calculations to data must be integrated over 
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b (fm) b (fm) b (fm) 

FIG. 14: (Color online) The suppression factor -RdPb ai y — —4 (left), (center) and 4 (right) as a function of b. The result is 
shown for J/ip (top) and T (bottom) in d+Pb relative to p + p collisions at the same energy, ys„„ = 6.2 TeV, and employ 
the EKS98 (solid), nDSg (dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. The horizontal lines 
show the impact-parameter integrated results. 



the width of the impact parameter bin which will average the impact parameter dependence of the shadowing over 
the bin width. Our calculations include a width of 0.2Ra for the impact parameter bins. 

In fact, studying Rcp in p + A and d-|-^ coUisions could provide a direct measure of shadowing if absorption 
is neghgible since higher-order corrections unrelated to shadowing cancel in the ratio 62]. As an example of an 
asymmetric system, Fig. [1^] presents Rcp{y) for d-|-Pb collisions with he = ^ and bp k, Ra^ . As expected, the 
resulting Rcp{y) are very similar to the impact-parameter averaged {d+ A) / {jp + p) ratios shown in Fig. 1^1 Since 
Rcp{y) with hp K, 2Ra are similar to those in Fig. [161 they are not shown. 

Figures flTl and ITSl show the values of Rcp for 6p « Ra and ^Ra relative to &c = in the four A-\- A systems 
studied for the J/i/' (Fig. [iTl) and the T (Fig. [HI). Since the change in i?AA(A^part) between 6c = and bp w Ra is 
small (see Fig. [TH]), these ratios are almost independent of rapidity and give Rcp close to unity. On the other hand, 
the weaker shadowing effect at hp « 2Ra produces a stronger rapidity dependence and a lower Rcp- Note that, 
as in Fig. [THl Rcp{y) for bp « 2Ra is similar to [{A + A)/{p + p)y>o in Fig. [TT] Thus, if no other medium effects 
are present, it is possible to trace the shadowing effect rather accurately by determining Rcp for sufficiently narrow 
centrality bins. 



III. SUMMARY 



We have provided a survey of the quarkonium total cross sections to next-to-leading order in the color evaporation 
model for all ^ -I- i? combinations and energies at the LHC. We have included initial-state shadowing, employing 
several parameterizations of the nuclear modifications of the parton densities, but assumed final-state absorption is 



^ We do not show Rcp for p + A collisions. 
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FIG. 15: (Color online) The suppression factor Raa at j/ = as a function of A'^part. The effect of shadowing on J ftp (a, 
upper 4 panels) and T (b, lower 4 panels) production is shown. Each set of panels displays the suppression factor for O+O at 
= 7 TeV (upper left), Ar+Ar at y/^^ = 6.3 TeV (upper right), Sn+Sn at y^s^ = 6.14 TeV (lower left) and Pb+Pb at 
y^s^ = 5.5 TeV (lower right). The calculations are with CTEQ6 and employ the EKS98 (solid), nDSg (dashed) and EPS09 
(solid curves with symbols) shadowing parameterizations. 



negligible. If the nuclear absorption of quarkonium production can indeed be ignored at LHC energies, it may be 
possible to use the different mass scales for J/ ip and T production to study the scale dependence of the gluon density 
in the nucleus as well as in the proton. There are considerable uncertainties in the predictions due to the incomplete 
knowledge of the nuclear gluon distribution. Indeed, at midrapidity, the range of the EPS09 {p + A)/{p + p) ratios 
differs by a factor of two. 

To illustrate the range of predictions for the different systems, we have calculated {p + A)/{p+p) and {d+A)/{p+p) 
ratios from the most naive (both systems at the same energy) to the most realistic (the p + p reference at 14 TeV and 
the rapidity shift oip + A interactions in the equal-speed frame). The most naive ratios are most straightforward for 
extracting the nuclear gluon distributions. It is still possible to use the most realistic ratios using a combination of 
experimental cuts on the rapidity distributions when Ayl^ ^ and modeling the appropriate x values for comparing 
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FIG. 16: (Color online) The central-to-peripheral ratios, Rcp, as a function of rapidity for bp ~ Ra relative to fe = for d+Pb 
collisions at ^/s^ = 6.2 TeV. The calculations are with CTEQ6 and employ the EKS98 (solid), nDSg (dashed) and EPS09 
(solid curves with symbols) shadowing parameterizations. 
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FIG. 17: (Color online) The central-to-peripheral ratios, Rcp{y), for b = iiA (a) and b = 2_Ra (b) relative to 6 = 0. The effect 
of shadowing on J/tp production is shown for O-l-O at ^/sj^ = 7 TeV (upper left), Ar-|-Ar at y^sj^ — 6.3 TeV (upper right), 
Sn-I-Sn at ^/sj^ = 6.14 TeV (lower left) and Pb-I-Pb at ^/sj^ = 5.5 TeV (lower right). The calculations are with CTEQ6 and 
employ the EKS98 (solid), nDSg (dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. 
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FIG. 18: (Color online) The central-to-peripheral ratios, Rcp{y), for b = Ra (a) and b = 2Ra (b) relative to b = 0. The effect 
of shadowing on T production is shown for 0+0 at y^s^^iv ~ TeV (upper left), Ar-|-Ar at ^Js^^ = 6.3 TeV (upper right), 
Sn+Sn at ^s^^ = 6.14 TeV (lower left) and Pb+Pb at ^s^^ = 5.5 TeV (lower right). The calculations are with CTEQ6 and 
employ the EKS98 (solid), nDSg (dashed) and EPS09 (solid curves with symbols) shadowing parameterizations. 



p + p collisions at -^s = 14 TeV with lower energy p -I- A collisions. As is clear from the RHIC analyses [l^, the A + A 
studies require a good understanding of the nuclear gluon distribution to extract hot and dense matter effects. Thus 
the more realistic d-l- A scenario shown in Fig. [10] would be preferable for determining nuclear effects on the parton 
densities both because of the relatively similar d+A and A + A center-of-mass energies as well as the smaller rapidity 
shift relative to p + A collisions in the equal-speed frame. 

To more cleanly extract the parton densities at LHC energies, it would be preferable to have e + p and e + A 
data at the appropriate x and fi^ range of the LHC data. (The HERA x range reaches to approximately the value 
appropriate for J/?/; production in 5.5 TeV/nuclcon collisions at midrapidity. Unfortunately, the /i^ probed at these 
X values is smaller than the J/tp mass scale.) So far, the nDIS data is not available at small enough x values and, 
simultaneously, large enough /i^ to be relevant for quarkonium production at high energies. While electron-proton 
collisions, as studied at HERA, would be useful for obtaining the baseline in p + p, they are not sufficient for defining 
the modification of the nuclear gluon distributions for p + A collisions, e + A studies are needed. 

The shadowing parameterizations used in our study exhibit a wide range of behavior for the nuclear gluon density 
at low X, outside the current range of the fits from fixed-target nDIS data at higher x and low . If nuclear data were 
available from high energy e + A collisions, the nuclear gluon densities could be more precisely pinned down by global 
analyses of the scale dependence of the nuclear structure functions. In hadroproduction, direct photon or open charm 
production, dominated by gluon-induced processes but without the additional complexities of nuclear absorption, 
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could be utilized to study the nuclear gluon density. Any new e + A data before an electron ring is available at the 
LHC will be at lower energies than previously available at HERA, reducing the potential overlap of the low x range 
between an electron- ion collider and the LHC. 

We note that, since we have assumed absorption is negligible at the LHC and include no other cold nuclear matter 
effect, the uncertainties on the ratios can be obtained from the EPS09 bands shown in the figures. However, if other 
effects are incorporated, a more extensive error analysis, including the uncertainties on these other effects, is necessary. 

Finally, we note that the central-to-peripheral ratio, Rcp, may be useful for extracting the shadowing effect at a 
given collision energy if the experimental resolution of the impact parameter bins is narrow enough. Indeed, Rcp 
measurements may be a superior method of studying asymmetric systems since very peripheral collisions are a good 
approximation to nucleon-nucleon collisions. This ratio is advantageous because it can be made at the same collision 
energy with a common rapidity shift. 
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